clc
clear
L = 10;
N = 128;
N2 = N / 2;
multi =  2 * pi * L; 
w = zeros(N,N);




for i=1:N
    for j = 1:N
       w(i,j) = sqrt( (multi * (j-N2))^2 + (multi * (i-N2))^2);  
    end
end
%w = fftshift(w);
omega = @(t) exp(-1i * w *t);

%x =  fftshift( full(createPsiK2D([2 2],[0.125 0.125],[0 0],L,N)));
x =   full(createPsiK2D([8 0],[0.05 0.05],[-4 5],L,N));
t = 0;
f = figure
while true
    if ishandle(f)
        t = t + 0.0001;
        next = x.*(omega(t));
        imagesc(abs(ifft2(next)))
        drawnow
    else;
        break;
    end;
end;
